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Abstract 

We use a simple picture based on the tt electron approximation to study the 
bandgap variation of carbon nanotubes with uniaxial and torsional strain. We 
find (i) that the magnitude of slope of bandgap versus strain has an almost 
universal behaviour that depends on the chiral angle, (ii) that the sign of 
slope depends on the value of (n — m) mod 3 and (iii) a novel change in sign 
of the slope of bandgap versus uniaxial strain arising from a change in the 
value of the quantum number corresponding to the minimum bandgap. Four 
orbital calculations are also presented to show that the tt orbital results are 
valid. 
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I. INTRODUCTION 



The mechanical and electronic properties of carbon nanotubes (CNT) have individually 
been studied in some detailEhi and the predicted dependence of bandgap on chiralityEHi has 
been observed.! The study of bandgap variation with mechanical deformation is important 
in view of the ability to manipulate individual nanotubesi. Additionally, they could form 
the basis for nanoscale sensors in a manner similar to experiments using Cqo molecules.^ 
Recent studies of bandgap change of zigzag and armchair tubes on mechanical strain have 
shown interesting behavior .Hill Refs. |S] and |H] studied the effect of uniaxial strain using a 
Green's function method based on the 7r electron approximation and a four orbital numerical 
method, respectively. Ref. |TT] predicted the opening of a bandgap in armchair tubes under 
torsion, using a method that wraps a massless two dimensional Dirac Hamiltonian on a 
curved surface. In this paper, we present a simple and unified picture of the band gap 
variation of chiral and achiral CNT with uniaxial and torsional strain. The method used 
is discussed in section |I[ The results obtained by using a single n orbital are discussed in 
section [111 A] and are compared to four orbital calculations in section n.rB[ The conclusions 
are presented in section [TV]. 



II. METHOD 



In the presence of a uniform uniaxial and torsional strain, a distorted graphene sheet 
continues to have two atoms per unit cell [Fig. 1]. It is convenient to represent the change 
in bond lengths using the chirality dependent coordinate system. The axes of the chirality 
dependent coordinate system corresponding to (n, m) CNT are the line joining the (0, 0) 
and (n, m) carbon atoms (c), and its perpendicular (i) [Fig. The fixed and chirality 

dependent coordinate system are related by, c = cosOx+sinOy and i = — sin9x+cos6y, where, 
sin{6) = \ IL ^ L and cos{6) = ^^^p- Ch = \/n 2 + m 2 + nm, is the circumference of the tube 
in units of the equilibrium lattice vector length, | The bond vectors are given 
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by, 

a n + m „ n — m~ ao m ~ 2n + m~ 

n = — c = 1 + dr 1 and r 2 = — c H -= 1 + dr 2 , (1) 

2 c ft 2V3 c h 2 c ft 2v/3 c h 

where, 5fi represents deviation from an undistorted sheet and r*3 = — (r*i + r 2 ). Within 
the context of continuum mechanics, application of a uniaxial or torsional strain causes the 
following change in the bond vectors of Fig. 1: 

r it — ► (1 + e t )r it and r ic — > (1 + e c )r ic (tensile) (2) 
-> r ic + tan(^f)r it (torsion), (3) 

where, i — 1,2,3 and r ip is the p-component of f\ (p = c,t). e t and e c represent the strain 
along t and c respectively, in the case of uniaxial strain. 7 is the shear strain. 
Using Eqs. (|])-(|3]), the lattice vectors of the distorted sheet are, 

-> r/-, A2n + m \^3m \f?>m~ 

ai = n - r 3 = a [(l + e c J- h tan{^)— Jc + a (l + e t )— 1 (4) 

^ r/ , 1 n + 2m . -V^n^ ,.. . v^3 n - ._. 
a 2 = 7*1 - r 2 = a [(l + e c )- tan(~f)— Jc - a (l + e t ) — 1 . (5) 

Z C/j Z C/j z c^ 

The corresponding unit cell area is |a x x a 2 | = ^(1 + e 4 )(l + e c )aQ. The real space unit cells 
correspond to r = j\di + j 2 a 2 , where ji and j 2 are integers. The ID unit cell length (T) is 
the shortest r t for which r c = 0. That is, the two lattice points, r = and r = jiSi + j 2 a 2 
have the same c coordinate. This corresponds to the following condition on ji and j 2 , 

(1 +e c )[ji(2n + m) + j 2 (n + 2m)] +tan(7)v / 3[jim-j 2 ^] =0 (6) 

and the ID unit cell length is, 

rp /, , v v 7 ^ (jim - j 2 n) 

T = a (l + e t )— (7) 

When only uniaxial strain is present (7 = 0), Eq. (0) corresponds to, j\{2n + m) + 
j 2 (n + 2m) = 0. The corresponding j\ and j 2 with smallest absolute values are j\ = 
(n + 2m) / gcd(2n + m,n + 2m) and j 2 = — (2n + m)/gcd(2n + m, n + 2m). gc<i refers to the 
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greatest common divisor. Using these values in Eq.(^), the ID unit cell length of an (n,m) 
tube is, 

T = (1 + e t )v3c^ao/ gcd(2n + m,n + 2m) . (8) 

In the absence of strain, Eq. (|8[) reduces to the result for undeformed nanotubes. In the 
presence of uniaxial strain, the unit cell length is equal to (1 + e t ) times the unstrained unit 
cell length. When only torsion is present, Eq. @ simplifies to, 

ji(2n + m) + ]2(n + 2m) + tan(7)v3(jim — j^n) = . (9) 

For arbitrary values of 7, n and m, this equation corresponds to a large T. For example, 
from Fig. 1 it is easy to see that under torsion, the unit cell of an armchair tube can be 
much larger than a$ depending on the value of 7. We will come back to this point at the 
end of section Q, where we discuss calculation of bandgap change due to torsion. 

We treat the nanotube within the approximation that it is a rolled up graphene sheet 
and assume a single it orbital per carbon atom. We calculate the band structure of the 
distorted sheet to be,0 

E(k) = {t\ + t 2 2 + t 2 3 + 2ht 2 cos[k ■ (fi - r 2 )] + 2t 2 t 3 cos[k ■ (f 2 - r 3 )] + 2t 3 ti cos[£ • (f 3 - , (10) 

where, k = k c c + k t t. The primary effects of change in bond vectors are to alter the hopping 
parameter between carbon atoms and lattice vectors. The hopping parameter is assumed to 
scale with bond length as,0 tj = to ( r o/ r i) 2 ) where to and r are the hopping parameter and 
bond length of an unstrained graphene sheet. The value of to is around 3eV. From Eqs. 
and (|j), the circumference of the distorted sheet is (1 + e c )c/ l a . The wave function of the 
CNT is quantized around the circumference and so k c is given by, 

fc c (l + e c )c h a = 2nq , (11) 

where, q is an integer. Eq. ( |T0"D can now be written as, 
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2 2 2 n + 2m a/3 n , . a/3 tanOy) n 

h{k t ) = {tj +t 2 + t 3 + 2*1*2 cos[7rg ^ q ^t a o — — j 

2n + m v^3^,/ \/3 £071(7) m 

+ 2*1*3 cos|7rg 2 1 K a o + 



2 c fc 1 + e c 4 

n — m y/3n + m,. v3 tanM n + m 1 
+ 2* 2 t3 cos[vrg — + ^ ^a + 7rg V w; — ])3 , (12) 



where, fcj = (1 + et)fc*. The bandgap of an (n, m) tube in presence of uniaxial (7 = 0) or 
torsional strain (e c = e t — 0) can be easily calculated from Eq. (Jl2|). In case of uniaxial 
strain, the limits of k t are given by — ^ < kf < ^, where T is the ID lattice vector length 
determined by Eq. The number of atoms in the ID unit cell does not change in the 
presence of uniaxial strain and so the range of q does not change from the undeformed case 
(q = 0, 1, 2, N c , where N c is the number of hexagons in the ID unit cell). 

In the case of torsion, the number of atoms in the ID unit cell and T can be large [Eq. 
(§)]. The corresponding span of k t is then small compared to the undeformed tube and the 
range of q is commensurate with the number of atoms in the ID unit cell. The eigen spectrum 
can however be obtained from Eq. (|I2"D by setting 7 = and spanning over the same values 
of q and k t as in the undeformed case. This is because the eigen spectrum depends only on 
the tight binding parameters (and not on the exact geometry) if the coordination number 
of the carbon atoms remains constant.0 

III. RESULTS AND DISCUSSION 

The results obtained using the method described in section |TJ are discussed in section 



I II A| . We then present the results from four orbital calculations with energy minimized 



structures in section f l [B 



A. 7r orbital 



We first consider the case of uniaxial strain. The bandgap is obtained by finding the 
minimum of E(kt), where the span of k t and q are discussed below Eq. (|TJ). The bandgap 
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change is largest for zigzag tubes and the magnitude of \dE g /da\ is approximately equal to 
3t . For armchair tubes, application of uniaxial strain does not cause a bandgap. We find 
that, (i) \dE g /da\ increases with increase in chiral angle [Fig. 2] and (ii) the sign of dE g /da 
follows the (n — m) mod 3 ruleS For example, the chiral angle of (6, 5) and (6, 4) tubes are 
close to that of armchair tubes. The slope of bandgap versus strain is correspondingly small 
and the sign of slope are opposite. For semiconducting zigzag tubes and armchair tubes, 
our results agree with Ref. |[ 

As uniaxial strain increases, there is an abrupt reversal in sign of dE g /da as illustrated 
for zigzag tubes in Fig. 4. This feature indicates a change in band index q corresponding to 
the bandgap and can be understood from the following expression that describes dependence 
of energy for various values of q at k t = [Eq. (|16D of appendix]: 



where, sgn(x) = [1 — 2cos(qir/n)]. The minimum value of E (q) = t |1 — 2cos(g7r/n)| is half 



for the band index q = q that satisfies n = 3qo ± 1. The second term can however change 
sign when q changes from qo to qo ± 1. As a result, a dramatic change in sign of dE g /da 
becomes possible if the magnitude of the second term is larger than change in the first term 
[Fig. 3] . The strain required to observe this effect decreases as the inverse radius of tube for 
large n. This is because the difference in energy of the qo and qo ± 1 bands become smaller 
with increase in radius. Fig. 3 demonstrates this point by comparing the (10,0) and (19,0) 
tubes. For the (19,0) tube, the change in slope occurs at around five percent strain. These 
strain values are accessible in bulk nanotube samples.0 The inset of Fig. 3 shows change 
in energy of the q=6 and q=7 bands for the (19,0) tube for three different values of strain. 
While the q=6 band shifts up in energy as strain increases, the q=7 band shifts down. Thus 
leading to the discussed change in sign of dE g /da. 

In case of torsional strain, the bandgap is obtained by finding the minimum of E(k t ) 
using Eq. (0), where the span of k t and q are discussed in the last paragraph of section ||. 




(13) 



of the bandgap of an unstrained tube. The first term of Eq. ( |T3"D takes the smallest value 
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The magnitude of \dE g /da\ is approximately equal to 3to for armchair tubes and this is in 
agreement with Ref. |TT[ For zigzag tubes, torsion causes only a small change in bandgap. 
The leading term in bandgap change depends on 7 only to second order. We find that, (i) 
\dE g /da\ decreases with increase in chiral angle and takes the smallest value for zigzag tubes 
[Fig. 4] and (ii) the sign of dE g /da follows the (n — m) mod 3 rule.0 

B. Four orbital 

To verify the simple picture presented, we have also performed four orbital calculations 
using the parametrization given in Ref. [18|. The change in bond lengths are computed using 
both continuum mechanics [Eqs. (0) and (|3|)] and structures that are energy minimized 
by Brenner potential!!! The energy minimization was performed with periodic boundary 
conditions. For the small values of strain considered, we find that the bandgap is not very 
sensitive to the two methods of obtaining the bond lengths. The results presented in Figs. 5 
and 6 correspond to the bond lengths obtained by energy minimization. For semiconducting 
tubes, the results of Figs. 5 and 6 agree with the 7r orbital results presented in Figs. 2 and 
4 respectively: The slope of dE g /da follows the (n — m) mod 3 rule and the magnitude of 
slope varies monotonically with chiral angle. The primary difference concerns non armchair 
tubes satisfying n — m = 3 * integer. This is not surprising because Ref. ^| has predicted 
such tubes to have a small bandgap due to curvature induced hybridization at zero strain. 
As a result, applying either tension or compression does not produce the "V" shaped curve 
of Fig. 2 with zero bandgap at zero strain. The difference is that the curves are shifted away 
from the origin as shown in Fig. 5. 

IV. CONCLUSIONS 

In conclusion, we present a simple picture to calculate bandgap versus strain of CNT 
with arbitrary chirality. We find that under uniaxial strain, \dE g /da\ of zigzag tubes is 3to 
independent of diameter, and continually decreases as the chirality changes to armchair, 
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when it takes the value zero. In contrast, we show that under torsional strain, \dE g /da\ 
of armchair tubes is 3to independent of diameter, and continually decreases as the chirality 
changes to zigzag, where is takes a small value. The sign of dE g /da follows the (n—m) mod 3 
rule in both cases.0 We also predict a change in the sign of dE g /da as a function of strain, 
corresponding to a change in the value of q that corresponds to the bandgap minimum. 
Comparison to four orbital calculations show that the main conclusions are unchanged. The 
primary difference involves nonarmchair tubes that satisfy n — m = 3 * integer. 

This work is supported by NASA contracts NAS2-14031 to Eloret (LY), NASA Ames 
Research Center and U.S. Department of Energy (JPL). 

V. APPENDIX 



Zigzag tubes under tension: Under uniaxial strain the band structure of (n, 0) is, 
E(k t ) = ±t 2 



1 jl / 4 *i\ t ( l' n \ rn , \ V^3 2ti 2 2 qn 

1 ± ( — ) cos( — ) cos((l + e t J— Mo) + (— ) cos ( — ) 

t2 n 2 C2 Tl 



t\ = t$ due to symmetry. The minimum of E(kt) occurs at kt = 0, 



E(0) = ±t 2 \l- — coa(^)\ . 

t 2 n 



To first order in 5ri Eqn. (15) is,0 



E(0) = E (q) - 2t, 



Sri 
r 



25r 2 qn 
1 - -=— cos — 
or i n 



sgn[x) 



1 - 2cos(^) 
n 



(14) 



(15) 



(16) 



where, sgn(x) = [1 — 2cos(g7r/n)] and E (q) = t |1 — 2cos(g7r/n)|. 
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Figure Captions: 

Fig. 1: The fixed (x,y) and chirality dependent (c,t) coordinates, ri, r 2 and r 3 correspond 
to bonds 1, 2 and 3 respectively. a\ and a 2 are the lattice vectors of the two dimensional 
sheet. 

Fig. 2: Bandgap versus tensile strain: For semiconducting tubes, the sign of slope of 
d(Bandgap) / d(Strain) depends only on the value of (n — m) mod 3. The magnitude of 
d(Bandgap) / d(Strain) is largest for zizag tubes and decreases with decrease in chiral angle. 
The magnitude is smallest for armchair tubes. The solid, dashed and dotted lines correspond 
to (n — m) mod 3 values of 1, -1 and respectively. The value of t is around 3eV. 

Fig. 3: The change in slope of the (10,0) and (19,0) tubes around 10% and 5% strain 
respectively is due to a change in the quantum number q that yields the minimum bandgap. 
Inset: E vs k of the q=7 (solid) and q=6 (dashed) bands as a function of strain for a (19,0) 
tube. Strains of 0%, 3% and 6% correspond to increasing thickness of the lines. 

Fig. 4: Bandgap versus torsional strain: For semiconducting tubes, the sign of slope of 
d(Bandgap) / d(Strain) depends only on the value of (n — m) mod 3. The magnitude of 
d(Bandgap) / d(Strain) is largest for armchair tubes and decreases with increase in chiral 
angle. The magnitude is smallest for zigzag tubes. The solid, dashed and dotted lines 
correspond to (n — m) mod 3 values of 1, -1 and respectively. 

Fig. 5: Same as Fig. 2 with the only difference that these are four orbital results. In the 
y-axis label, t=2.66eV. 

Fig. 6: Same as Fig. 4 with the only difference that these are four orbital results. In the 
y-axis label, t=2.66eV. 



1-11 



FIGURES 




Fig. 1 / Yang 



1-12 




STRAIN (%) 



Fig. 2 / Yang 



1-13 





ng 4 / Yang SHEAR (in degrees) 



1-14 



0.6 



XT 0.5 

O 

</> 

.tS 0.4 



^ 0.3 
Q_ 
< 
O 

Q 



0.2 



< 

CQ 0.1 



0.0 



(8,0) 


































* 


» 




" ■ — - 


- — « 






















~ — < 


► 




































- 


" ■ 4 


— _^ 
















(7 5 


) 

1 







- - 


1 


i 






• 




X 











































































































































a 

h 




























i 

i 


k 

k 


h0) 



















1 






* 


k 




°>V 

5,5) 


— i 


k 




Si 






...-ik 

"... — Jk 
- — Ik 



-4 -3 -2-10 1 

STRAIN (%) 



Fig. 5 / Yang 



0.6 



^ 0.5 

O 

(/) 

.tr 0.4 



CL 
< 
O 

Q 



(9,1) 



03 w 



0.2 



< 

CD 0.1 



0.0 



(Si) 



(5,5) 



-4.0 -3.0 -2.0 -1.0 0.0 1.0 2.0 

Fig 6 / Yang SHEAR (in degrees) 



3.0 4.0 



1-15 



